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ABSTRACT 



Context. At least four high-frequency quasiperiodic oscillations (QPOs) at frequencies 41 Hz, 67 Hz, 113 Hz, and 167 Hz were re- 
ported in a binary system GRS 1915+105 hosting near-extreme Kerr black hole with a dimensionless spin a > 0.98. 
Aims. We attempt to explain all four observed frequencies by an extension of the standard resonant model of epicyclic oscillations. 
Methods. We use the idea of oscillations induced by the hump of the orbital velocity profile (related to locally non-rotating frames- 
LNRF) in discs orbiting near-extreme Kerr black holes, which are characterized by a "humpy frequency" v h , that could excite the 
radial and vertical epicyclic oscillations with frequencies v r , v v . Due to non-linear resonant phenomena, the combinational frequen- 
cies are allowed as well. 

Results. Assuming mass M = 14.8 M G and spin a = 0.9998 for the GRS 1915+105 Kerr black hole, the model predicts frequencies 
Vh = 41 Hz, v T = 67 Hz, Vh + v r = 108 Hz, and v v - v r = 170 Hz corresponding quite well to the observed ones. 
Conclusions. For black-hole parameters being in good agreement with those given observationally, the forced resonant phenom- 
ena in non-linear oscillations, excited by the "hump-induced" oscillations in a Keplerian disc, can explain high-frequency QPOs in 
near-extreme Kerr black-hole binary system GRS 1915+105 within the range of observational errors. 

Key words. Accretion, accretion discs - Black hole physics - Methods: data analysis 
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1. Introduction 

Detailed analysis of the variable X-ray black-hole binary system 
(microquasar) GRS 1915+105 reveals high-frequency QPOs ap- 
pearing at f our frequencies, namely v\ - (41 + l )Hz, V2 = 
', (67 ± l)Hz dMorgan et alJll997t IStrohmaveill200lh . and v 3 = 
(113 + 5)Hz, v 4 = (167 + 5)Hz lRemillardl2004 . In this range of 
its errors, both pairs are close to the frequency ratio 3:2 suggest- 
ing the possible existence of resonant phenomena in the system. 
Observations of oscillations with these frequencies have differ- 
ent q ualities, but in all four cases the data are quite convincing; 
see (IMcClintock & Remillardl 120041 IRemillard & McClintockl 
120061) . 



Several models have been developed to explain the kHz QPO 
frequencies, and it is usually preferred that these oscillations are 
related to the orbital motion near the inner edge of an accretion 
disc. In particular, two ideas base d on the strong-gravity prop- 
erties have been proposed. While IStella & Vietril (U998l 1 19991) 
introduced the "Relativistic Precession Model" considering that 
the kHz QPOs directly manifest the modes of a slightly per- 
turbed (and therefore epicyclic) relati vistic motion of blobs in 
the inn er parts of the accretion disc, iKluzniak & Abramowiczl 
J2001 propose models based on non-linear oscillations of an 
accretion disc that assume resonant interaction between orbital 
and/or epicyclic modes. In a different context, the possibility 
of resonant coupling between the epicyclic modes of motion 



in the Kerr spacetime w as also mentioned in the early work of 
lAliev&Galtsovl (119811) . 

In the case of near-extreme Kerr black holes, it was sug- 
gested that the epicyclic oscillations in the disc could be ex- 
cited by resonance s with the so-called "hum p-in duced" oscilla- 
tions, s ee pap ers of lAschenbachl d2004l 12006k and lStuchlfk et al.l 
d2^04Tl200% . This idea was proposed so as to extend standard 
orbital (resonant) models meant to explain high-frequency QPOs 
observed in black-hole sources. 

Recently, careful and detailed analysis of the spectral con- 
tinuum from GRS 1915 + 105 has put a strong limit on the black- 
hole spinQ namely 0.98 < a < 1 jMcClintock et alj|2006l) . indi- 
cating the presence of near-extreme Kerr black hole whose mass 
has been restricted observationally t o M = (14.0 ± 4.4) M , 
see (IMcClintock & Remillardl l2004t IRemillard & McClintockl 
120061) . Therefore, the microquasar GRS 1915+105 seems to be 
an appropriate candidate to test the extended resonant model 
with hump-induced oscillationsQ 

The idea of hump-induced oscillations and their possible res- 
onant coupling with the epicyclic ones is briefly discussed in 
Sect. [2] The related resonant model, assuming the excitation of 
epicyclic oscillations by the hump-induced oscillations through 
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1 Units c = G = M = 1 (M is the total mass of the Kerr black 
hole) and the Boyer-Lindquist (B-L) coordinates (t, r, 9, <p) are used 
hereafter. 

2 However, Middleton et al. ( 2006) refer to a substantially lower, in- 
termediate value of black-hole spin, a ~ 0.7, to which the model of 
hump-induced oscillations cannot be applied. 
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non-linear resonant phenomena, is applied to GRS 1915 + 105 in 
Sect. [3] concluding remarks are presented in Sect. [4] 

2. Hump-induced and epicyclic oscillations in 
Keplerian discs and possible resonant coupling 

In order to describe the local processes in an accretion disc, 
it is necessary to choose a local observer (characterized by its 
reference frame). In general relativity there is no preferred ob- 
server. On the other hand, if we want to study processes related 
to the orbital motion of matter in the disc, it is reasonable to 
choose the observers with zero angular momentum, so-called 
ZAMOs, as their reference frames do not rotate with respect to 
the spacetime, and thus ZAMOs should reveal local kinematic 
properties of the disc in the clearest way. (In rotating-stationary, 
axisymetric-spacetimes, they are dragged along with the space- 
time.) In the Kerr spacetime, ZA MOs are represented b y locally 
non-rotating frames (LNRF); see iBardeen et alj d!972l) . Notice 
that in the Schwarzschild spacetime, LNRF correspond to the 
sta tic observer frames. 

lAschenbachl d2004l) finds that for near-extreme Kerr black 
holes with the spin a > 0.9953, the test-particle orbital veloc- 
ity < V^ > related to LNRF reveals a hump in the equatorial plane 
(0 = n/2). This non-monotonicity is located in a small region 
inside the ergosphere of the black-hole spacetime close to, but 
above, the marginally stable orbitQ Therefore, it can be relevant 
for thin accretion discs around near-extreme Kerr black holes, 
as the inner edge of the disc can extend down to the innermost 
stable circular orbit (ISCO). 

Moreover. [Stuchlfk et alj d2005l) shows that for a > 0.99979 
the similarly humpy behavior of the orbital velocity in LNRF 
also takes place for the non-geodesic motion of test perfect fluid 
in marginally stable barotropic tori characterized by the uni- 
form distribution of the specific angular momentum, £(r, 6) = 
-Uip/Ut = const, where the motion of fluid elements is given 
by the 4-velocity field V = (U'(r, 0), 0, 0, lP(r, 0)). Outside 
the equatorial plane, the non-monotonic behavior of T^' in 
marginally stable tori is represented by the topology change of 
the cylindrical equivelocity surfaces in the region of the hump, 
because the toroidal equivelocity surfaces centered around the 
circle corresponding to the local mi nimum of 'V^ in the equa- 
torial plane exist for a > 0.99979 dStuchlfket al.ll2005l) . This 
suggests a generati on of possible instabil ities in radial and ver- 
tical directions; see IStuchlfk et al.l d2004l) . In the following, we 
restrict our attention to the case of Keplerian discs. 

Heuristic connection between the positive part of the veloc- 
ity gradient, d^V^'/dr, and the excitation of epicyclic oscilla- 
tions in Keplerian discs was suggested by lAschenbachl (12004 
120061) . who defined the characteristic frequency of oscillations, 
induced by the humpy profile of <yW, by the maximum posi- 
tive slope of the orbital velocity in terms of the coordinate ra- 
dius, v cr it = (<yy (y) /c) r) max . This co o rdinat e-dependent defini- 
tion was corrected in IStuchlik et al. (2007), where the proper 
radial distance dr = yfg^dr rather than the coordinate dis- 
tance dr was used to define the characteristic (critical) frequency 
•/. = (d'V^ I <9r) max . Such a locally defined critical frequency 
was further related to a stationary observer at infinity, obtaining 
the so-called "humpy frequency" 
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= yj-(gtt + lugty + w 2 g w ),. = ,. h v' clit 

3 We stress that the Aschenbach effect is frame-dependent, as it is 
related to LNRF, but recall the arguments for relevance of the LNRF 
point of view at the beginning of the section. 
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Fig. 1. Spin-dependence of the humpy frequency Vh and the 
humpy radius r^. For completeness, the B-L radius of the in- 
nermost stable circular geodesic (ISCO) is plotted. 
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where g MV are the metric coefficients of the Kerr geometry and 
w — —ghflgifip is the angular velocity of the LNRF; see, e.g., 
IBardeen et alj d!972l) ; Ah = r^ - 2rh + a 2 . The analytic formula 
is given for the equatorial plane (0 = n/2). The B-L radius 
where the positive gradient of the velocity profile in terms of the 
proper radial distance reaches its maximum, so-called "humpy 
radius", is given by the condition 



dr \ df 







(2) 



leading to the equation 

3a 7 (r + 2) + a 6 VK21r 2 + 18r - 4) - a 5 r(33r 2 + lOr + 20) 
+a 4 r V?(45r 3 - 62r 2 - 68r + 16) - aV(83r 2 - 122r - 60) 
+a 2 r 4 VK27r 2 - 130r + 136) - 9ar 5 (7r 2 - 26r + 24) 
+r 7 VK3r - 2) = 0, (3) 

which must be solved numerically. The spin dependence of the 
humpy radius and the related humpy frequency is illustrated in 
Fig. Q] The humpy radius falls monotonically with increas- 
ing spin a, while the humpy frequency Vh has a maximum for 
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a - 0.9998, where v h ( m ax) = 607 (Mq/M) Hz, and it tends to 
VtC^D =588 (Mq/M) Hz. 

When particles following a Keplerian circular orbit are 
perturbed, they begin to follow, in the first approximation, 
an epicyclic motion around the equilibrium Keplerian or- 
bit, generally characterized by the fre quencies of the radial 
and vertical epicyclic oscillations v r , v v dAliev & Galtsov|[T981t 
lNowak&LeM1998l) : 



4d 



3aV 2 ), 



6r l + 8ar~ 3/2 
v; = v 2 = v 2 (1 - 4ar- y2 + 3a V 2 ), 
where vk is the Keplerian orbital frequency 



1 



v K = 



2jr(r y2 + a) ' 



(4) 
(5) 

(6) 



The ratios of the humpy frequency and th e epicyclic freq u encies 
at the humpy radius were determined in IStuchlfk et al.l d2007l) 
revealing almost spin-independent asymptotic behavior for a — > 
1 represented closely by the ratios of integer numbers, v v : v r : 
Vh ~ 1 1 : 3 : 2, which imply a possibility of resonant phenomena 
bet ween the hump-ind uced and epicyclic oscillations predicted 
by lAschenbachl d2004l) . The ratios of the epicyclic frequencies 
and the humpy frequency are given in the dependence on the 
black-hole spin in Fig. [2] 
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Fig. 2. Spin dependence of frequency ratios including the ra- 
dial (v r ) and vertical (v v ) epicyclic frequencies, and the humpy 
frequency (vh) evaluated at the same radius where the 
humpy frequency is defined. The range of the spin relevant for 
GRS 1915+105 Kerr black hole is shaded. For the mean value 
a = 0.9998, the frequency ratios are close to the ratios of in- 
teger numbers, suggesting a possibility of resonances between 
hump-induced and epicyclic oscillations in GRS 1915+105. 



3. Application of the hump-induced resonance 
model to high-frequency QPOs in GRS 1915+105 

Primarily concentrating on the lower pair of frequencies, we as- 
sume that the lowest frequency is directly the humpy frequency, 

v h = v 1 = (41±l)Hz, (7) 

while the second lowest frequency corresponds directly to the 
radial epicyclic frequency at the same radius r^, 

v 1 - = v 2 = (67 + l)Hz. (8) 

These frequencies are close to a 3:2 ratio, therefore the forced 
non-linear resonance can be relevant in such a situation. The 
ratio of v r /vh = (1.63 + 0.06) gives the black hole spin a = 
(0.9998 + 0.0001) (the uncertainty of the spin is implied by un- 
certainties of the lower pair of frequencies being ~ 1 Hz); see 
Fig. [2] Notice that this spin corresponds to the maximal pos- 
sible value of the humpy frequency Vh(max) (Fig. [T). Since the 
humpy frequency is 1/M-scaled, the absolute value of im- 
plies the black hole mass M = (14.8 + 0.4) Mq. The correspond- 
ing humpy radius is = 1.29 ^ JJ' (Fig. [TJ. At such a radius, 
the vertical epicyclic frequency of a particle orbiting the Kerr 
black hole with the mass and spin inferred above reaches the 
value v v = (0.23 + 0.01) kHz. Then the upper pair of observed 
frequencies can be explained, within the range of observational 
errors ±5 Hz, by combinational frequencies at the humpy radius 
r-n in the following way: 

v 3 ~ (v r + v h ) = (108 + 2) Hz (9) 
v 4 ~ (v v - v r ) = (0.17 ± 0.01)kHz. (10) 

4. Conclusions 

The idea of epicyclic oscillations induced by the LNRF-velocity 
hump in the region where the positive part of the velocity gra- 
dient reaches its maximum is able to address all four high- 
frequency QPOs observed in the X-ray source GRS 1915+105. 



The model implies a near-extreme spin of the central black hole 
{a ~ 0.9998), which agrees well with results from the spec- 
tral continuum fits, and the black-hole mass M ~ 14.8 Mq be- 
ing well inside the interval given by other observational meth- 
ods. Note that the orbital resonance model of Kluzniak & 
Abramowicz, assuming the parametric resonance between the 
vertical and radial epicyclic oscillations in frequency ratio 3:2 
represented by the upper pair of obser ved frequencies, a lso gives 
the spin a > 0.99 but for M 18M Q dTorok et alj|2005l) . On the 
other hand, the "Relativistic Precessio n Model" gives a s ubstan- 
tially lower value for the spin: a ~ 0.3 dStella et al.lll999l) . 

In the presented model, we assume that all four observed 
frequencies arise due to forced non-linear oscillations of the 
Keplerian disc at the same radius r^, excited by the hump- 
induced oscillations characterized by the humpy frequency Vh. 
The black-hole parameters a, M are fixed by the requirement 
that the lower pair of observed frequencies is identified with the 
humpy frequency and the radial epicyclic frequency, v\ = Vh, 
V2 = v r . Assuming non-linear resonant phenomena enabling 
the existence of combinational frequencies and the possibility 
of observing them, the upper pair of observed frequencies can 
be explained as the combinational ones of the humpy frequency 
and both epicyclic frequencies, V3 ~ (v r + Vh), v\ ~ (v v - v r ). 
Moreover, both frequency ratios v r : Vh, and (v v - v r ) : (v r + Vh) 
are close to 3:2 ratio (Fig. |2j, in which the resonant phenomena 
can be strong enough. On the other hand, as 4vh = (164 + 4) Hz, 
which is also close to the uppermost frequency, there is an- 
other possibility of explaining V4 through a sub-harmonic reso- 
nance forced by the humpy oscillations as well. Finally, note that 
IStrohmaverl d2001l) also reports another relatively weak QPO at 
frequency of (56 + 2) Hz. If this is the case (which, according 
to our knowledge, has not been confirmed by other observations 
yet), it could be related to the second harmonic of the combina- 
tional frequencjQ (v r - Vh) = (26 ± 2) Hz. 



4 Combinational frequency (v r - Vh) corresponds to the same order of 
nonlinearity as (v r + Vh). 

Note added in the manuscript: After the paper was accepted we ob- 
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Generally, other harmonics and combinational frequencies 
may occur in a non-linear oscillating system corresponding to 
higher approximations, when the equation of motion describing 
the non-linear oscillations is solv ed by the method of succe ssive 
approximations. The statement bv lLandau & Lifshitzld 19761) that 
"As the degree of approximation increases, however, the strength 
of the resonances, and the widths of the frequency ranges in 
which they occur, decrease so rapidly that in practice only the 
resonances at frequencies v ss pvo/q with small p and q can be 
observed" can explain why a QPO near the frequency 237 Hz, 
corresponding to the vertical epicyclic frequency v v at the same 
radius as the previously mentioned humpy and radial epicyclic 
frequencies Vh, v r , is not directly observed, despite the commen- 
surability of these frequencies represented by the frequency ra- 
tios v v : Vh ~ 6 : 1 and v v : v r ~ 7 : 2 (Fig.|2]i. 
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ta ined an information t hat a weak QPO at frequency 27 Hz is referenced 
iri lBelloni etai] feOOlh . 
5 p, q are integers. 



